The potential of the ground state of NaRb 
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The X 1 E+ state of NaRb was studied by Fourier transform spectroscopy. An accurate potential 
energy curve was derived from more than 8800 transitions in isotopomers 23 Na 85 Rb and 23 Na 87 Rb. 
This potential reproduces the experimental observations within their uncertainties of 0.003 cm -1 to 
0.007 cm -1 . The outer classical turning point of the last observed energy level (v" = 76, J" = 27) 
lies at ta 12.4 A, leading to a energy of 4.5 cm -1 below the ground state asymptote. 

PACS numbers: 31.50.Bc, 33.20.Kf, 33.20.Vq, 33.50.Dq 



I. INTRODUCTION 

The heteronuclear alkali dimers attract interest of both 
experimental and theoretical researchers involved in colli- 
sion dynamics, photoassociative spectroscopy, laser cool- 
ing and trapping of alkali atoms [J, Q, EJ l| ■ Special inter- 
est is put on the study of the ground states and especially 
near the atomic asymptote, since the precise knowledge 
of the long-range interactions between two different types 
of alkali atoms is necessary for understanding and real- 
ization of such cold collision processes as sympathetic 
cooling, formation of two species BEC (TBEC) and ul- 
tracold heteronuclear molecules. 

Although the NaRb molecule is one of the promising 
candidates for TBEC Q, the experimental information 
on the ground singlet state is still limited. 

Accurate term values of the X 1 E + state were deter- 
mined by the Doppler-free polarization spectroscopy 
They cover a range of low vibrational quantum numbers 
v" G [0, 6] with limited sets of rotational quantum num- 
bers (e.g. J" e [1,80] for v" = and J" G [18,36] for 
v" = 6). Later Kasahara et al. @ extended the range of 
v" up to v" = 30, but only for J" = 10, 12. On the other 
hand the recent experiment by the Riga group @ col- 
lected about 300 transition frequencies to energy levels 
with a much wider spectrum of vibrational and rotational 
quantum numbers (v" G [15,76], J" G [12,64]) but with 
relatively moderate accuracy ( 0.1 cm^ 1 ). A hybrid po- 
tential for the NaRb ground state based on extrapolation 
between regions determined from the early experimental 
data and the theoretical dispersion coefficients from 
Ref. |8( was constructed by Zemke et al. |9| giving an im- 
proved estimate on the dissociation energy of the ground 
state - D e = 5030.75 cm" 1 . 

On the theoretical side two papers have recently ap- 
peared: by Korek et al. 0] an d by Zaitsevskii et al. 
[llj j . In both of them theoretical potential energy curves 
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FIG. 1: Selected low singlet and triplet electronic states in 
NaRb according to Ref. |lC||. 



for the ground and several excited states are given. Re- 
sults from Ref. 10] on selected low electronic states are 
shown also in Figure ^ which will help in following the 
electronic assignment of the new observations, reported 
below. 

The goal of the present investigation is to study ex- 
tensively the ground electronic state X X E + of the NaRb 
molecule. We chose the Fourier transform spectroscopy 
which is able to provide abundant and precise spectro- 
scopic data, needed for construction of an accurate po- 
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tential energy curve. In addition, this study is considered 
as the necessary background for further investigations of 
this molecule at long internuclear distances required for 
modeling cold collisions at lmK and below. 



II. EXPERIMENT 

The NaRb molecules were produced in a single section 
heat-pipe oven similar to that described in Ref. 0] by 
heating 5 g of Na (purity 99.95 %) and 10 g of Rb (purity 
99.75 %, natural isotopic composition) from Alfa Aesar. 
The oven was operated at temperatures between 560 K 
and 600 K and typically with 2 mbar of Ar as buffer 
gas. At these conditions apart from the atomic vapors 
all three types of molecules are formed - Na2, Rb2 and 
NaRb. This mixture was illuminated by three different 
laser sources: a single mode, frequency doubled Nd:YAG 
laser, Ar + ion laser and Rhodamine 6G dye laser. Their 
wavelengths fall in a spectral window with weak absorp- 
tion of Rb2- Setting the working temperatures as low 
as possible we reduced the density of the Na2 molecules 
reaching a relatively intense NaRb spectra free of Na2 
emission. The heat pipe oven was in operation more than 
60 hours without refilling and at the end of the experi- 
ment was still in good working condition. 

The fluorescence from the oven was collected in the 
direction opposite to the propagation of the laser beam 
and then focused into the input aperture of Bruker IFS 
120 HR Fourier transform spectrometer. The signal was 
detected with a Hamamatsu R928 photomultiplier tube 
or a silicon diode. The scanning path of the spectrometer 
was set to reach a typical resolution of 0.0115 cm -1 -j- 
0.03 cm" 1 and the number of scans for averaging varied 
between 10 and 20. In order to avoid the illumination of 
the detector by the He-Ne laser, used for calibration of 
the spectrometer, a notch filter with 8 nm FWHM was 
introduced. For better signal-to-noise ratio some spectra 
were recorded by limiting the desired spectral window 
with colored glass or interference filters. 

The single mode, frequency doubled Nd:YAG laser 
with a typical output power of 70 mW at 532.1 nm ex- 
cited the C 1 E+ <- X 1 S+ and D 1 !! <- X 1 E+ systems 
of NaRb, compare Fig. 1. The frequency was varied 
between 18787.25 cm" 1 and 18788.44 cm -1 and spectra 
were recorded at frequencies which excite strong fluores- 
cence. 

The Ar + ion laser was operated both in single (typi- 
cal power 100-^500 mW) and multi mode (typical power 
0.5-^3 W) regimes. The 514.5 nm, 501.7 nm, 496.5 nm, 
488.0 nm and 476.5 nm Ar + ion laser lines induced flu- 
orescence mainly due to the D 1 !! «— X 1 S + system of 
NaRb. The Na 2 B X II U -> X : S+ band was also observed, 
especially exciting with the bluer Ar + laser lines, but its 
intensity was reduced by decreasing the working temper- 
ature down to 560 K. Along with the D <— X system, 
the 514.5 nm line excites also transitions in the C 1 S + <— 
X 1 E + system of NaRb Q- The special form of poten- 
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FIG. 2: The vibrational progression up to v" = 76 in Na 85 Rb 
excited by a single mode Ar + laser. Lines not indicated are 
also assigned. 



tial energy curve of the C 1 S + state 01 results in favor- 
able Franck-Condon factors for transitions from some C 
state levels both to low and high lying energy levels of 
the ground state. Such C state levels are excited by the 
514.5 nm line and the green florescence of the D <— X 
system is accompanied by a red fluorescence due to C <— 
X transitions, which reach v" — 76 of the ground state 
(Fig.©. 

In order to bring the frequency of the single mode 
Ar + laser exactly in resonance with desired transitions 
we did the following. As a first step we recorded the 
LIF for a given frequency simultaneously with the Fourier 
spectrometer and a GCA/McPherson Instruments scan- 
ning monochromator (lm focal length). The high reso- 
lution Fourier spectrum helped us to assign unambigu- 
ously the record of the monochromator. Then we set the 
monochromator on a strong and pronounced line of the 
progression of interest and we tuned the Ar + laser fre- 
quency until the maximum signal on the monochromator 
was reached. The frequency stability of the laser with 
a temperature stabilized intracavity etalon was sufficient 
in order to perform 10 to 20 scans with the Fourier spec- 
trometer. 

The Rhodamine 6G laser excites the B 1 !! <— X 1 E + 
transition and also weak transitions in the C 1 S + <— 
X 1 S + system and was used for two reasons. The first 
was to enrich the information on the ground state lev- 
els with intermediate vibrational quantum numbers be- 
cause of the gap between the levels from the D x n <— 
X 1 !]^ system (mainly with low v") and those from the 
C 1 S + <— X 1 !] 4 " system (mainly with high v"). Secondly, 
we wanted to excite levels of the B state with significant 
triplet admixture due to perturbations from the neigh- 
boring b 3 LT and c 3 S + states 0, 0. Indeed, scanning 
the frequency of the dye laser from 16729 cm" 1 to 16965 
cm" 1 we encountered a large number of excitations where 
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FIG. 3: The [B^, b 3 n, c 3 E+]-X 1 E+ and [B 1 !!, b 3 n, c 3 E+]- 
a 3 E + bands of NaRb excited by a single mode Rhodamine 6G 
laser. The insert shows the continuum due to the bound- free 
emission to the repulsive wall of the a 3 E + state. 



a second fluorescence band around 12000 cm -1 appeared 
along with the B <— X system (see Fig. |3J). The analysis 
of this band (a report is in preparation) confirmed that 
we have observed transitions to the ground triplet state 
a 3 E+. 

All transitions excited in Na 85 Rb and Na 87 Rb by the 
single mode Nd:YAG and Rhodamine 6G lasers and by 
the Ar + laser in single- and multi mode, assigned and 
used in this analysis are listed in Table 1 of the supple- 
mentary material (EPAPS No ). 



III. ANALYSIS 

The assignment of the recorded spectra was simplified 
by the published potential energy curve in Ref . Q • After 
we identified the strongest progressions, which were al- 
ready observed in previous studies 0, HH Hj, Ha] , a new 
potential was fitted. This preliminary potential was fur- 
ther improved filling additionally collected experimental 
data; this sequential procedure allows continuous check- 
ing of the assignment, especially for high J" and in cases 
of large gaps in v". The total data set (Fig. consists 
of more than 6150 transitions in Na 85 Rb and 2650 in 
Na 87 Rb given in Table 2 of the supplementary material 
(EPAPS No ). 

For fitting the experimental data we used the approach 
adopted in Ref. [T3| - In order to exclude any fitting 
parameters concerning the excited states, we fitted the 
ground state potential directly to the observed differences 
between ground state levels. In our case of 8800 transi- 
tions we selected about 43300 differences for representing 
the different cases of vibrational intervals (for details see 
also 01 )• The uncertainties of the differences were esti- 



FIG. 4: The range of vibrational and rotational quantum 
numbers of the observed ground state energy levels in Na 85 Rb 
and Na 87 Rb 

mated taking into account the instrumental accuracy and 
resolution, and the signal-to-noise ratio for both transi- 
tions forming the difference. The uncertainty of strong 
lines was set to 0.002 0.003 cm -1 increasing to 0.007 
cm -1 for weaker ones with SNR w 2. Differences which 
were measured several times were averaged with weight- 
ing factors determined from their uncertainties. 

As already discussed in [l^, for our experimental set 
up a fluorescence progression excited by a single mode 
laser will be shifted in frequency if the laser is not tuned 
to the center of the Doppler profile. This problem is par- 
tially overcome by fitting differences instead of transition 
frequencies. Nevertheless there can be still a residual 
shift of the differences due to the frequency dependence 
of the Doppler shift. For example if the laser frequency 
around 20000 cm -1 is detuned from the line center by 
one Doppler width (about 0.03 cm -1 for NaRb at 600 
K) the difference between two transitions v\ and vi sep- 
arated by v\ — V2= 4000 cm" 1 will be shifted from the 
"true" one by 0.006 cm" 1 . This shift is comparable with 
the typically expected uncertainty (w 0.004 cm" 1 ) but 
it cannot be determined in our experiment. Therefore, 
initially we did not correct the uncertainties of the large 
differences (a possible solution could be to increase the 
uncertainty of the difference v\ — V2 proportionally to 
the ratio (y\ — v?)lv\). Large differences however were 
checked carefully during the fitting procedure and finally 
no need for such correction was found. 



IV. CONSTRUCTION OF PEC 

The ground singlet states of both isotopomers of NaRb 
are described in the adiabatic approximation with a sin- 
gle potential energy curve. The potential was constructed 
as a set of points {Ri,U(Ri)} (see Ref. [16]) connected 
with cubic spline function. The pointwise representation 
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TABLE I: Pointwise representation of the potential energy 
curve for the X X E+ state of NaRb. 
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FIG. 5: The derived standard deviation (in logarithmic scale) 
of the points of the pointwise potential curve for three dif- 
ferent values of the singularity parameter Q (for details see 
Ref. E2). 



of the potential, however, is not convenient for long inter- 
nuclear distances. On the other hand, in order to ensure 
proper boundary conditions for solving the Schrodinger 
equation for high vibrational energy levels, distances up 
to R w 20 A are needed. Therefore for R > i? out we 
adopted the usual dispersion form: 



U(R) = D e - C 6 R~ e - C S R- 8 - C W R- 



(1) 



with coefficients Cg , Cs taken from Ref. |18l llPJ | . The 
connecting point R ou t and the parameters D e and C\q 
were varied in order to ensure a smooth connection with 
the pointwise potential as follows. Initially, the PEC was 
constructed in a pointwise form up to 16 A. After each 
fitting iteration D e and C\o were adjusted to fit the shape 
of the pointwise potential between 13.0 A and 16.0 A to 
the analytic formula The crossing point of the point- 
wise and the long range curves was taken as i? ut • Since 
the experimental data are almost insensitive to details of 
the potential shape in this region, we found this proce- 
dure sufficient and did not fit the long range parameters 
directly to the experimental data. Once the final form 
of the PEC was achieved, we found that it is possible 
to extend the validity of the long range expression down 
to approximately 11.8 A. At the same time the quality 
of the derived potential did not change if the analytic 
formula was fitted to the shape of the pointwise po- 
tential between 11.8 A and 13.24 A. A discussion on the 
advant age of using this matching procedure can be found 
in Ref. |3- 

The final set of potential parameters consists of 51 
points. In order to calculate the value of the potential 
for R < R ou t a natural cubic spline 0J through all 51 
points and for R > R ou t the long range expression (JTJ 
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should be used with parameters listed in Table [I] This 
final PEC gives a standard deviation of a = 0.0031 cm -1 
and a normalized standard deviation of a = 0.70 showing 
the internal consistency of the data set and the quality 
of the derived potential. 

In order to find the region where the PEC is unambigu- 
ously characterized by the experimental data we analyzed 
the uncertainty of the fitting parameters as described in 
detail in Ref. 0, |2l|. By minimizing the merit func- 
tion x 2 using the Singular Value Decomposition (SVD) 
method flU it is possible to determine the parame- 
ters on which \ 2 and consequently the quality of the fit 
depend only weakly. The fitting parameters in our case 
are the values of the potential itself in a preset grid of 
points. In Fig. [3] the uncertainties of the fitted potential 
for intermediate internuclear distances are presented for 
three different values of the singularity parameter £ (see 
|2l|). Although the outer classical turning point of the 
last observed energy level (v" = 76, J" = 27) is 12.4 A 
we see that the shape of the potential energy curve is un- 
ambiguously fixed (within the experimental uncertainty) 
approximately only between 3.0 A and 9.8 A since the 
uncertainties there almost do not depend on £. There- 
fore, although the fitted potential energy curve describes 
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all the experimental data up to v" = 76, there is no rigor- 
ous way to estimate the uncertainty of its shape beyond 
9.8 A (see also 12]). The result of Fig.Elmight indicate 
that the correlation between the left (R < 3A) and right 
(R > 9.8 A) branch of the potential becomes significant. 
In this respect the reported potential in Table [I] for these 
regions is only one possibility which describes our obser- 
vations within the error limits. From Table [I] we derive 
some parameters which will be useful when applying the 
potential for other calculations and spectroscopic studies. 
Their values are not rounded in order to keep consistency 
with the data from Table HJ 

• equilibrium distance of the potential, 
R e = 3.643415 A, 

• position of the level with v" = and J" = 0, 

i?oo = 53.3117 cm -1 with respect to potential min- 
imum, 

• the dissociation energy of the ground state, calcu- 
lated with respect to Eqo-, 

.0^=4977.5363 cm" 1 . 

In Fig. the potential energy curves published in 
Ref. 0, 13 are compared with the pointwise potential 
(PW) of this study. The reason for the difference be- 
tween the PW potential and the Modified Lennard-Jones 
(MLJ) potential around R e is probably that for the 
construction of the later potential the rotational quan- 
tum numbers for levels with v" < 30 were limited only 
to the narrow interval from 6] (J" = 10, 12). With such 
data set the rotational constant B v can be determined 
only with large uncertainty which corresponds to shifts 
of the whole potential along the internuclear axes which 
is clearly indicated in Fig.[B] As already discussed in Sec- 
tion [I] the hybrid potential from Ref. Q was constructed 
using the whole range of available rotational quantum 
numbers from the literature (|5j, |6j) and therefore the 
agreement with the present study around R e is much 
better. The deviations between the PW and the hybrid 
potential reaching 20 cm -1 for intermediate internuclear 
distances are caused by the extrapolation between the 
theoretical long range part of the potential and the ex- 
perimental short range part. In the same region the ex- 
perimentally determined shape of the MLJ potential is 
in a much better agreement with the present study. 

Since the primary data from [5j were not available to 
us, we were able to check our potential only against the 
data published in [6J. Forming differences between tran- 
sition frequencies exactly as for our LIF data we found an 
agreement with the differences predicted by the present 
PEC from Table U with a standard deviation of 0.0026 
cm" 1 , well within the expected experimental error of 100 
MHz. Only two transitions ((12,11)-(17,10) and (12,11)- 
(19,10)) were excluded from this comparison since differ- 
ences with them showed systematically large deviations 
reaching 0.01 cm" 1 . 
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FIG. 6: Differences between the pointwise representation of 
the NaRb ground state potential derived in this study (PW) 
and previous potentials energy curves from Ref. 9] (Hybrid), 
Ref. (MLJ). The energy at the equilibrium point R e is set 
to zero for all potentials and used as a reference point. 

In TablelTTla set of Dunham coefficients for 23 Na 85 Rb is 
listed which describes the whole set of the measured dif- 
ferences for both isotopomers with a deviation of 0.0030 
cm" 1 and a normalized standard deviation of 0.84. For 
calculating the eigenvalues for 23 Na 87 Rb the usual scal- 
ing rules were used. The distribution of the indexes of the 
Dunham coefficients indicates that for description of such 
a large set of vibrational and rotational quantum num- 
bers the parameters of the Dunahm expansion lose their 
original meaning of "spectroscopic constants". The ex- 
tremely small values of some coefficients (~ 10~ 40 cm" 1 
!) indicate that they should be used with caution due to 
roundoff errors. At last, the extrapolation properties of 
such set of Dunham coefficients are doubtful. Therefore, 
we prefer to give the full description of the experimen- 
tal data by a potential energy curve, which is obviously 
a much better physical model, and to present Dunham 
coefficients only for convenience since there are still ap- 
plications where their use is easier and faster. 

V. RESULTS AND DISCUSSION 

In this paper a spectroscopic study of the NaRb ground 
state is presented. A variety of laser sources was used to 
excite several band systems in both isotopomers of the 
NaRb molecule and the induced fluorescence was ana- 
lyzed with a high resolution Fourier transform spectrom- 
eter. More than 8800 transitions were identified resulting 
in a set of more than 4000 ground state energy levels with 
wide range of vibrational and rotational quantum num- 
bers. An accurate potential energy curve was fitted to 
the experimental data reproducing differences between 
ground state energy levels with a = 0.0031 cm" 1 and 
a = 0.70. 

In order to check the consistency of the derived poten- 
tial energy curve we performed additional fits in several 
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TABLE II: Dunham coefficients for 23 Na 85 Rb in cm 1 , derived from the experimental data with v" < 76 and J" < 200. 



i 


k 






1 





106.8544964 




2 





-0.37984464 




3 





-7.703864 


xlO" 4 


5 





-5.390011 


xlO" 7 


7 





6.5551619 


xHT 10 


9 





-1.9725192 


xlO" 12 


10 





9.9663375 


xlO" 14 


11 





-2.5123 


xlO" 15 


12 





3.7418 


xlO" 17 


13 





-3.33869 


xlO" 19 


14 





1.652 


xlO" 21 


15 





-3.48 


xlO" 24 





1 


7.019021 


xlO" 2 


1 


1 


-2.943389 


xlO" 4 


2 


1 


-1.699605 


xlO" 6 


4 


1 


-1.474183 


xlO" 9 


6 


1 


5.732233 


xlO" 12 


7 


1 


-3.68622 


xlO" 13 


8 


1 


1.031766 


xlO" 14 


9 


1 


-1.2551 


xlO" 16 


11 


1 


1.22142 


xlO" 20 


13 


1 


-1.6097 


xlO" 24 


14 


1 


9.446 


xlO" 27 





2 


-1.20357 


xlO" 7 


1 


2 


-7.98199 


xlO" 10 


2 


2 


-4.21837 


xlO" 12 


7 


2 


-1.2888117 


xlO" 17 


8 


2 


1.760426295 


xlO" 18 


9 


2 


-1.17937534 


xlO" 19 


10 


2 


4.906556 


xlO" 21 


11 


2 


-1.3413739 


xlO" 22 


12 


2 


2.3980927 


xlO" 24 


13 


2 


-2.684354 


xlO" 26 


14 


2 


1.70 


xlO" 28 


15 


2 


-4.6349 


xlO" 31 





3 


1.4803 


xlO" 13 


1 


3 


1.9481 


xlO" 15 


3 


3 


-9.3543 


xlO" 18 


8 


3 


7.87129 


xlO" 24 


9 


3 


-6.789833 


xlO" 25 


10 


3 


2.13759 


xlO" 26 


11 


3 


-2.61852 


xlO -28 


13 


3 


1.974 


xlO" 32 


15 


3 


-1.0266 


xlO" 36 


1 


4 


-1.953 


xlO" 20 


2 


4 


-1.906 


xlO" 21 


5 


4 


1.11596 


xlO" 24 


6 


4 


-7.344 


xlO" 26 


8 


4 


7.4791 


xlO" 29 


10 


4 


-6.261 


xlO" 32 


11 


4 


8.982 


xlO" 34 


14 


4 


-6.3065 


xlO" 40 


7 


5 


-1.3005 


xlO" 33 



7 



different ways. Instead of fitting differences, we applied 
the commonly used approach of fitting directly transition 
frequencies, getting the term energies of the excited state 
levels as free parameters. In this case two representations 
of the potential were tried: a pointwise and an analytic 
form described in detail in Ref. p^ . Both representa- 
tions gave the same quality of the fit compared to the 
pointwise potential determined from differences and the 
potentials were practically identical in the studied range. 
Thus no additional parameters sets are given here in or- 
der to avoid the confusion of the reader what potential 
form should be preferred. 

The classical turning point of the last observed energy 
level (v" = 76, J" = 27) is around 12.4 A. Although 
this point lies beyond the Le Roy radius for the NaRb 
3 2 S+5 2 S asymptote (i?LeRoy = 11-2 A), the analysis of 
the PEC uncertainties indicates that with the present 
body of experimental data it is unsafe to determine Cq 
and the other dispersion coefficients. For this reason 
the experimental PEC was smoothly connected to a long 
range potential formed by dispersion terms taken from 
the literature 0, . For better description of the near 
asymptotic region of the PEC transition frequencies to 
weakly bound energy levels are needed (see for exam- 
ple [20]). Some improvement might be expected already 
from the combined analysis of the X 1 E + and a 3 E + states 
at long intcrnuclear distances which is in progress. 

Weiss et al. Q applied the MLJ potential from and 
extended it for R > 11 A by the long range parameters 
from 0, and the exchange interaction taken from [24[ 
to determine the singlet and triplet scattering lengths for 
3 2 S Na + 5 2 S Rb. Using the newly derived potential for 



the inner part we constructed the full potential with the 
same long range parameters as in |4| and derived a singlet 
scattering length of +305ao (do ~ 0.52918 A, the Bohr 
radius) for 23 Na 85 Rb and +103a for 23 Na 87 Rb. These 
results do not agree with the ones from Ref. Q: +167ao 
(+50a /-30a ) for 23 Na 85 Rb and 55a (+3a /-3a ) for 
23 Na 87 Rb. We used exactly the same approach for the 
long range as in Q to see clearly the influence of the 
inner part of the potential. By the new values we will not 
claim that we have determined more realistic scattering 
lengths. Rather we want to point out that the existing 
data do not allow to specify the scattering length in an 
useful uncertainty interval. Note also that slight changes 
in the exchange energy still allow to alter the sign of the 
scattering length in 23 Na 85 Rb. Thus energy levels closer 
to the asymptote (< 4.5 cm -1 ) must be studied. 
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